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Abstract 



In this paper, we investigate the instability of one-dimensionally stable periodic traveling wave 
£f} . solutions of the generalized Korteweg-de Vries equation to long wavelength transverse perturbations 

in the generalized Zakharov-Kuznetsov equation in two space dimensions. By deriving appropriate 
C^) \ asymptotic expansions of the periodic Evans function, we derive an index which yields sufficient 

■ conditions for transverse instabilities to occur. This index is geometric in nature, and applies to any 

periodic traveling wave profile under some minor smoothness assumptions on the nonlinearity. We 
also describe the analogous theory for periodic traveling waves of the generalized Benjamin-Bona- 
|/~) ■ Mahony equation to long wavelength transverse perturbations in the gBBM-Zakharov-Kuznetsov 

' equation. 

ON ' 
O 

1 Introduction 

' It is well known that the Korteweg-de Vries (KdV) equation 

SH ; 

u t = u xxx + uu x , (i,t)el 2 

arises as an model for one-dimensional long wavelength surface waves propagating in weakly nonlinear 
dispersive media, as well as the evolution of weakly nonlinear ion acoustic waves in plasmas [23] . There 
are several weakly two-dimensional variations on the KdV equation. Of particular importance is the 
Zakharov-Kuznetsov (ZK) equation [26] 

u t = (u xx + u yy ) x + uu x , 

which arises when in the study of geophysical fluid dynamics in anisotropic settings and ionic-acoustic 
waves in magnetized plasmas. 

"The author would like to thank Jared C. Bronski for suggesting this problem and for his many useful conversations 
during the preparation and revision of this manuscript. Also, the author gratefully acknowledges support from a National 
Science Foundation Postdoctoral Fellowship under grant DMS-0902192. 
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In many applications, it is useful to introduce a parameter a in the nonlinearity as follows: 

u t — u xxx + OiUU x . 

Notice that by scaling one can always assume that a has a fixed sign. In the context of plasmas, there 
are situations in which the physical parameters in the problem (such as temperature and density) force 
a to vanish. In neighborhood of such parameter values, it is necessary to consider higher terms in the 
nonlinearity and one usually resorts to using a modified KdV (mKdV) equation 

Ut = Uxxx + (3u 2 u x 

to describe such ion acoustic waves. Here, the sign of /3 can not be forced to be definite by scaling 
considerations and the two signs of (3 correspond to different physical phenomenon. For example, the 
focusing mKdV equation (corresponding to f3 > 0) can be derived as a model for the evolution of ion 
acoustic perturbations with a negative ion component, while the defocusing mKdV (corresponding to 
f3 < 0) can be derived as a model for the evolution of ion acoustic perturbations in a plasma with two 
negative ion components of different temperature. Moreover, the structure of the exponentially decaying 
solutions are radically different: for (3 > the solitary waves are homoclinic to a constant state, while for 
j3 < the solitary waves correspond to shock like solutions which are heteroclinic to opposite constant 
states. As in the case of the KdV equation, the modified ZK (mZK) equation 

U t = {U XX + Uyy) x + (iU 2 U X 

occur naturally as weakly two-dimensional variations of the mKdV. 

In order to encompass as many physical applications as possible, the focus of our study will be on the 
generalized KdV (gKdV) equation 

u t = u xxx + f(u) x (1) 
along with the generalized ZK (gZK) equation 

Ut = (uxx + u yy ) x + f(u) x (2) 

For suitable nonlinearities (satisfying mild convexity and smoothness assumptions), the gKdV equation 
admits asymptotically constant solutions known as the solitary waves. It seems natural to question the 
stability of such solutions of the gKdV as one-dimensional solutions in the gZK equations. It is clear the 
no sense of stability can occur unless the solitary wave is itself stable as a solution of the corresponding 
y-independent equations. It is well known (see [5], [7], and [20]) that such planar solitary waves are 
one-dimensionally (orbitally) stable to localized perturbations as solutions of the gKdV if 

d f 

— / u 2 (x + ct; c)dx 
oc J R 

is positive, and are (exponentially) unstable to such perturbations if it is negative: here u is the profile 
of the solitary wave and c is the wave speed. Physically, this says that such a solution is stable precisely 
when its momentum is an increasing function of the wave speed. In the case of power-law nonlinearities 
f(u) = u p+1 , p > 1, this is equivalent with stability if 1 < p < 4 and instability if p > 4. 

The transverse instability of one-dimensionally stable solutions of the gKdV in the gZK equation has 
been well studied in the solitary wave context. In the case of the ZK equation, Bridges [9] derived a 
geometric condition for the long wavelength transverse instability of solitary solutions of ([I]) by using 
the multi-symplectic structure of the ZK. There, the author derived an index which could be viewed as 
the Jacobian of a particular map and whose sign determined the transverse stability of the underlying 
wave. Moreover, Bridges was able to evaluate this Jacobian explicitly for power-law nonlinearities f(u) = 
u p+1 /(p + 1) by using the fact that all solitary wave solutions in this case can be expressed in terms of 
hyperbolic functions. 
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If one rather considers the periodic solutions of ((T|), then the stability theory is much less understood. 
This reflects the fact that the spectrum of the corresponding linearized operators is purely continuous, 
and hence it seems more difficult for nonlinear periodic waves of the gKdV to be stable than their solitary 
wave counterparts. Moreover, the periodic waves of the gKdV in general have a much more rich structure 
than the solitary waves: even in the case of power-law nonlinearities it is not possible to write down a 
general elementary representative for all periodic solutions of the gKdV, which stands in contrast to the 
solitary wave theory. It should be noted, however, that the periodic solutions are perhaps more physically 
relevant as they are used to model wave trains: such patterns are prevalent in a variety of applications. 
Up to this point, there seems to be two types of stability results in the periodic case: nonlinear stability 
to perturbations with periods which are integer multiples of the underlying waves period (see [1], [2], [8], 
[12], and [19]), and spectral stability to localized perturbations (see [10], [14], and [16]). In general, the 
methods employed in the nonlinear stability analysis parallels that of the energy functional framework 
from the solitary wave theory [15]. However, such methods do not extend in a straight forward manner if 
one considers localized perturbations. Instead, the techniques in this case mainly center around studying 
the spectrum of the corresponding linearized operators and applying regular perturbation theory: see 
for example [16]. Recently [10], there has been progress with analyzing the linearized spectrum in a 
neighborhood of the origin in the spectral plane by using periodic Evans function methods introduced by 
Gardner [13]. Such unstable spectrum corresponds to instabilities to slow modulations of the underlying 
wave, i.e. to long wavelength perturbations. The authors in [10] essentially conduct a rigorous modulation 
theory calculation to derive asymptotic expansions of the Evans function near the origin. 

In this paper, we will extend the methods in [10] in order to analyze the spectral stability of y- 
independent solutions of the gZK to long wavelength transverse perturbations. In particular, we will 
conduct a rigorous modulation theory calculation. To this end we will first show the corresponding lin- 
earized spectral problem has a periodic eigenvalue at the origin of multiplicity greater than one. Using 
the integrability of the ODE defining the traveling wave profiles of © and analyzing the corresponding 
periodic Evans functions, we are able to analytically determine the leading order behavior of the periodic 
eigenvalues bifurcating from the origin in the wavenumber of the transverse perturbation. This infor- 
mation yields sufficient conditions for the spectrum of the perturbed spectral problems to have unstable 
spectrum in a neighborhood of the origin. These conditions are geometric in nature and can be expressed 
as the Jacobian of a particular map from the traveling wave parameters to the period and mass of the 
underlying wave. The positivity of this index implies modulational instability of the underlying periodic 
wave to long wavelength transverse perturbations. By using the results of [10], [11], and [17] we are able 
to evaluate the sign of this Jacobian in several cases involving power-law nonlinearities: in particular, we 
are able to make conclusive instability statements for the ZK and focusing mZK equations. 

Also of interest in this paper is the analogous periodic transverse instability theory for the generalized 
Benjamin-Bona-Mahony (gBBM) equation 



When f(u) = u 2 , equation (|3|) is precisely the Benjamin-Bona-Mahony equation, also known as the 
regularized long wave equation, which arises as an alternate model to the gKdV equation as a description 
of gravity water waves in the long wave regime (see [6] and [21]). As with the gKdV equation, a natural 
weakly two-dimensional generalization of ^ is the gBBM-ZK equation [25] 



Recently, the stability of periodic traveling waves of the gBBM equation to both periodic and localized 
perturbations has been conducted [18] using the Evans function techniques of [10]. In particular, it was 
found that the stability theory closely resembles that of the gKdV equation. As to be expected then, if one 
considers the transverse instability of such a solution of ([3]) to long wavelength transverse perturbations 
in the gBBM-ZK equation the resulting theory parallels that of the gKdV equation described above. 
Therefore, we will only outline the corresponding transverse instability theory in this case. 

The outline of this paper is as follows. In section 2 we recall from [10] and [17] the basic properties 
of the periodic traveling waves solutions of (JTJ) . In particular, we will discuss a parametrization of such 



tit - u xxt + u x + f(u) x = 0. 



(3) 



U t - (U xt + Uyy) x +U X + f(u) x = 0. 



(4) 
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solutions which is useful in our calculations. We then carry out our transverse instability analysis for the 
gKdV equation in section 3. In section 4, we outline the analogous theory for the gBBM equation and 
end with some discussion and closing remarks in section 5. 

2 Properties of the Periodic Traveling Waves of gKdV 

In this section, we recall the basic properties of the periodic traveling wave solutions of (JTJ) . A traveling 
wave solution of the gKdV with positive wave speed c > is a solution of the traveling wave ordinary 
differential equation 

U XXX + f(u) x ~ CU X = 0, (5) 

i.e. they are solutions of ([T]) which are stationary in the moving coordinate frame defined by x + ct. 
Clearly, such solutions are reducible to quadrature and satisfy 

u xx + f(u) - cu = a, (6) 
\u\ + F(u) - C -v? - an = E, (7) 

where a and E are real constants of integration, and F satisfies F 1 = /, F(Q) = 0. In order to ensure the 
existence of periodic orbits of f5J, we must require that the effective potential 

V(u; a, c) = F(u) — —u 2 — au 

has a non-degenerate local minimum. Notice this places a restriction on the admissible class of nonlin- 
earities as well as the allowable parameter regime for a given nonlinearity. This motivates the following 
definition. 

Definition 1. We define Q C M 3 to be the open set consisting of all triples (a,E,c) such that (|5|) has at 
least one periodic solution. 

When (a, E, c) G f2, we use the notation u(x; a, E, c) to denote a particular periodic solutiorQ- Given 
an (a, E, c) G O, we will always assume that the roots u± of E = V(u; a, c) which correspond to the 
absolute max and min, respectively, of our periodic wave are simple. It follows that u± are C 1 functions 
of a, E, and c on fl and, with out loss of generality, u(0) = u_. 

Remark 1. Taking into account the translation invariance of ([T]) 7 it follows that for each (a,E,c) € f2 
we can construct a one-parameter family of periodic traveling wave solutions of (jlj : namely 

itj(x, t) = u(x + ct + £;a,E,c) 

where £ £ M. Thus, the periodic traveling waves of (JXJ) constitute a four dimensional manifold of solutions. 
However, this added constant of integration does not play an important role in our theory and can be 
modded out. 



Given an (a, E, c) S f2, notice that the period of the corresponding solution can be represented as 

r u + du 



T = T(a,E,c) :=2 / 

J U- 



y/2 (E-V(u;a,c)) 



(8) 



The above integrand can be regularized at the square root branch points it_, u + by a standard procedure 
(see [10] or [17]). In particular, we can differentiate the above relation with respect to the parameters 



1 In the case where more than one such solution exists for a particular (a, E, c), we can distinguish them by their initial 
values. 
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a, E, and c within the parameter regime Q. Similarly the mass, momentum, and Hamiltonian of the 
traveling wave are given by the first and second moments of this density, i.e. 

[ T C u + 11 du 
M(a, E, c) — / u(x) dx = 2 (g) 



y/2 (E-V(u;a,c)) 



r 



P(a,E,c) = I u 2 (x) dx = 2\ U - (10) 

y/2{E- V(u;a,c)) 

H(a,E,c) = I ( ZS. - F(u) ) M dx = 2 P - ~ & ' c) — S du (11) 



^2 (E-V(u;a,c)) 

These integrals can be regularized as above, and represent conserved quantities of the gKdV flow re- 
stricted to the manifold of periodic traveling wave solutions. In particular one can differentiate the above 
expressions with respect to the parameters (a, E, c). 

Remark 2. Notice that in the derivation of the gKdV [3], the solution u can represent either the hor- 
izontal velocity of a wave profile, or the density of the wave. Thus, the functional M can properly be 
interpreted as a "mass" since it is the integral of the density over space. Similarly, the functional P can 
be interpreted as a "momentum" since it is the integral of the density times velocity over space. 

Throughout this paper, a large role will be played by the gradients of the above conserved quantities. 
Thankfully, the underlying classical dynamics of the |(SJ) provides many useful relations between these 
gradients. To see this, note that the classical action (in the sense of action angle variables) is given by 

u x du= u 2 x dx = 2 / \J2(E - V(u;a,c)) du. (12) 

JO Ju- 

The derivative of this map as a function K : f2 — > K is given explicitly by 

V a , E ,cK(a,E,c) = (^M(a,E,c),T(a,E,c),^P(a,E,c 

where V a ,E,c = (j§^)~§e>~§^)- This is easily verified by differentiating (fl2|) and recalling that E — 
V(u±:a,c) = 0. Using the fact that T, M, P, and H are C 1 functions of parameters (a,E,c) on 
the domain f2, we see that 

E V a ,E,c{T) + a V a.E.c 

see [10] for details. Thus, although the subsequent theory is developed most naturally in terms of the 
quantities T, M, and P, it is possible to restate our results in terms of M, P, and H so long as E ^ 0: 
this is desirable since these have a natural interpretation as conserved quantities of the partial differential 
equation (fT]). 

We now discuss the aforementioned parametrization of the periodic solutions of (|5|) in more detail. A 
major technical necessity throughout this paper, especially in the analysis of the gZK equation, is that 
the constants of motion for the PDE flow defined by (fT]) provide (at least locally) a good parametrization 
for the periodic traveling wave solutions. In particular, we assume for a given (a, E, c) € that the 
conserved quantities (H, M, P) are good local coordinates for the periodic traveling waves near (a, E, c) 
in the sense that the map (a,E,c) — > (H(a,E,c),M(a,E,c),P(a,E,c)) has a unique C 1 inverse in a 
neighborhood of (a, E, c). Adopting the notation 

r, hX _ d(f,g,h) 
d(x,y,z) 

for 3x3 Jacobians, we see that this is possible exactly when {H, M, P} a ,E,c is non-zero, which is equivalent 
to {T, M, P} a.E.c being non-zero if E ^ by (fT3"j) . In particular, it was shown in [10] that the vanishing 
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of {T, M, P} a ,E,c forces a change in the structure the generalized periodic null-space of the linearized 
operator: see Lemma [T] below for more details. This non-degeneracy condition has been shown to 
generically hold in several cases of power-law nonlinearities in [11]. In the case of the gZK equation, the 
sign of this Jacobian will play an important role in determining of the corresponding periodic traveling 
wave solution of the gKdV is a one-dimensionally stable solution of the gZK equation. When considering 
the transverse instability of such a solution, another important Jacobian which arises is {T, M} at E, 
where this 2x2 Jacobian is defined analogous to the above 3x3 case. Notice that the non- vanishing of this 
Jacobian implies the period and mass can be used to parameterize nearby periodic waves with a fixed 
wave speed. It is precisely the sign of this Jacobian that determines the stability of such a solution to 
long wavelength transverse instabilities in both the gZK and gKP equations. 

It should be noted that the geometric quantities {T,M} a ^E and {T, M, P} a ,E,c also arise naturally 
when considering the stability of a periodic traveling wave solution of the gKdV equation perturbations 
whose period is an integer multiple of that of the underlying wave (see [10], [11], and [17]). Since, as 
mentioned in the introduction, the long wavelength transverse instabilities occur from the bifurcation of 
periodic eigenvalues from the origin into the unstable space, it is not surprising that these quantities play 
a large role in the transverse instability analysis. 

3 Transverse Stability Analysis 

In this section, we begin our transverse stability analysis for one-dimensionally stable periodic traveling 
wave solutions of JT]). As mentioned before, a T-periodic traveling wave solution u(x,t) — u(x + ct) of 
([1]) corresponds to a y-independent T-periodic traveling wave solution solution of traveling gZK equation 

U t = (u xx +u yy ) x + f(u) x - cu x . (14) 

We now consider a small localized perturbation of u of the form 

tfj(x, y, t) = u(x) + ev(x, y, t) + 0{e 2 ) 

where v(x, •,*) € L°°(M) for each (x,t) G R 2 and v(-,y,t) <E L 2 (R) for each (y,t) € M 2 , and |e| < 1 is a 
small parameter. Forcing i/> to solve (|14l) yields a hierarchy of consistency conditions: the O(e ) equation 
holds since u is a solution of (fT4"|) . and the 0(e) equation reads 

-vt = d x (£[u] - dy) v 

where C[u] — —d 2 — f(u) + c. As this linearized equation is autonomous in time and in the spatial variable 
y, we may seek separated solutions of the forrrJl 

v(x,y,t) =e-» t+iky v(x) 

where fi € C, and k G M is the transverse wavenumber of the perturbation. This leads to the (ODE) 
spectral problem 

d x (C[u\ +k 2 )v = Aiw (15) 

considered on the real Hilbert space L 2 (M.). Notice that we consider C[u] as a closed, self-adjoint second 
order differential operator with periodic coefficients with dense domain in L 2 (M.) . 

Our general strategy in analyzing ()15|) is as follows. As the spectral problem involves differential 
operators with periodic coefficients, standard results in Floquet theory imply that the L 2 spectrum is 
purely continuous, and consists entirely of L°° eigenvalue^. We define the monodromy matrix M(/x, k) 

2 Alternatively, you could take the Laplace transform in time and the Fourier transform in the spatial variable y. 

3 Notice this shows one of the main difficulties in considering the spectral stability of periodic solutions as opposed to the 
asymptotically constant case. In the latter, the continuous spectrum can in general be characterized via a Weyl sequence 
argument and the Fourier transform and shown to not contribute to an instability. Thus, one must only study the L 2 point 
spectrum of the linear operators. In the periodic context, however, the spectrum is entirely essential and hence any spectral 
instability comes from the continuous spectrum. 
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to be the solution of the corresponding first order system 

Y x = H(x,n,k)Y, Y(0)=I, (16) 

where I is the three by three identity matrix. Following [13] we can then characterize the L°° eigenvalues 
as the roots of the periodic Evans function, which we define as 

D(fi, k, A) := det ( M(n, fc) - AI) . 

The complex constant A referred to as the Floquet multiplier and is related to the class of admissible 
perturbations considered: for example, A = 1 corresponds to periodic perturbations with the same period 
as the underlying wave. It can be shown that, for a fixed transverse wavenumber k, \x is in the spectrum 
of the corresponding linearized operator if and only if there exists a A € C with |A| = 1 such that 
D([i, k, A) = 0, i.e. if and only if the corresponding monodromy operator has an eigenvalue on the unit 
circle. As a consequence, the L 2 spectrum of the corresponding spectral problems, for a fixed k, can be 
parameterized by by the discrete zero-set of the functions D((a, k, A) with |A| = 1, i.e. it is equal to the 
set 

Moreover, it will be shown that the natural symmetries present in the gKP equations forces D(Q, 0, 1) = 0, 
and hence the origin in the spectral plane is a periodic eigenvalue for each of the spectral problems 
(|15p . Our goal is to understand the nature of the spectrum in a neighborhood of fi = when the 
transverse wavenumber k is very small, and hence we must study the zero set of the function D(fj,, k, 1) 
in a neighborhood of the point (fi,k,X) = (0,0,1). This will allow us to determine how these periodic 
eigenvalues bifurcate from the origin in the transverse wavenumber and will lead us to a sufficient condition 
for instability. 

We begin our analysis by studying transverse instability within the gZK equation. Throughout, we 
will work with a fixed periodic traveling wave solution u(x) a, E, c) of |T]) with (a, E, c) € Q and consider 
the transverse instability of u as a solution of (I14p . By linearizing (11411 around u and applying the 
appropriate transforms in time and the transverse direction y leads to the spectral problem 

d x [C[u] + k 2 ) v = fxv 

considered on L 2 (R). We begin by determining when a periodic traveling wave solution of the gKdV 
equation is a one-dimensionally stable solution of the gZK equation, i.e. determining necessary condi- 
tions which imply the L 2 (M) spectrum of the unperturbed linearized operator is confined to the 
imaginary axiqj. 

First, we recall the necessary results of [10] where the spectrum of the operator c^Z^m] was exten- 
sively analyzed in a neighborhood of /i = 0. The periodic traveling wave profile it(-; a, E, c) satisfies the 
traveling wave ODE and moreover satisfies the identities ^ and 0. In particular, it follows from 
differentiating (|6|) that the functions u x , u a , and ue are linearly independent and satisfy 

C[u]u x = 0, C[u]u a = 1, C[u]u E = 0. (17) 

It follows that we can use these three functions to form a basis to construct the monodromy matrix 
M(0, 0) in this case. In particular, it can be shown that the matrix M(0,0) has 1 as an eigenvalue 
of algebraic multiplicity three and geometric multiplicity two. Using this along with the Hamiltonian 
structure of the unperturbed linear operator c^/^u] and a bit of perturbation theory, we are able to prove 
the following lemma. 

Lemma 1. The function D(fi, k, 1) satisfies the following properties. 



4 It follows that by one-dimensionally stable, we mean that it is a y-independent solution which is stable to y-independent 
perturbations. 
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1. D(fi,k, 1) is an odd function of \x and an even function of fc. 

2. The limit lim^^oo sign {D{n, 0, 1)) exists and is negative. 

3. The asymptotic relation 

D((x, 0, 1) = -1{T, M, P} a , E , c fi 3 + 0(| M | 4 ) 

holds in a neighborhood of fi = 0. 

In particular, the corresponding periodic traveling wave solution of (JTJ) is one-dimensionally unstable as 
a solution of © if {T, M, P} a ,E,c < 0. 

Proof. The second and third parts have been proven in [10] using perturbation theory and asymptotic 
analysis. The fact that -D(/i, fc, 1) is odd in jj, was also proved in [10] using the fact that the Hamiltonian 
structure of the operator d x C[u] implies that M(/x, 0) ~ M(— ^,0) _1 . Similarly, L>(/i, fc,l) is even 
in fc since the perturbed spectral problem d x (C[u] + fc 2 ) v = fiv is invariant under the transformation 
k^-k. □ 

We are now interested in the case where the transverse wavenumber k is small. By Lemma [TJ fi = 
is a periodic eigenvalue of the linear operator d x C[u] of multiplicity three. We now wish to see how these 
periodic eigenvalues bifurcate from the origin in the transverse wavenumber k. To do so, we determine 
the dominant balance of the equation D(fi, k, 1) = for /i and k small. This is the content of the following 
lemma. 

Lemma 2. The equation /c,l) = has the following normal form in a neighborhood of (/i, k) = (0,0): 



-^-{T, M, P} a , E . c + fik 2 {T, M} a<E J u\dx + 0(H 4 + k 



4 ^-0. 



Proof. By standard arguments using a Newton diagram (see [4] or [24]), Lemma [T] implies the dominant 
balance of the equation D(/i, k, 1) = Q near (//, k) = (0, 0) must take the form 

-fj{T, M, P} a , E , c + ^-D„ kk {Q, 0, 1) + 0(H 4 + fc 4 ) = 0. 

Thus, it is enough to prove D^ki®, 0, 1) = 2{T, M} a ^E Jq u x dx. 
To begin, we write (fT5]) as a first order system of form (fTE|) with 

/ 10 

U(x, fi>,k)=\ 1 

V -fi-f"(u)u x c + k 2 -f'(u) 

From (JT7J), we can use the functions u x , u a , and ue to construct a corresponding matrix solution in the 
case (/i, fc) = (0, 0). To this end, we define ~W(x, fx, fc) to be the matrix solution which satisfies 

Ux Ua U E 

W(x, 0,0)= ( u xx u ax u Ex 

l&xxx ^axx W'Exx 

and fix the initial condition W(0,/i, fc) = W(0,0,0) for all (p,k) e C x M. Defining 6W(fi,k) 

I T 

W(x, fj,,k)\ Q , it follows that 

£>(/*,Jfe,l) = -det(5W(/i,A)) 
since det ( W(0, 0, 0)) = — 1. Thus, it suffices to analyze the matrix <JW(/x,fc). 
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Next, we compute 8 W(0, 0) with the plan of treating 8 W(^i, k) as a small perturbation of this constant 
matrix. Notice that by modeling out the translation invariance, we can set 

u(0) = u- = u(T) 
u x (0) = = u x (T) 
u xx {0) = -V'(u-) = u xx (T). 

Using these equations along with the chain rule and ([I]) , a straightforward calculation shows that 

/ 

£W(0,0)= V'{uJ)T a V'{u-)T E 
\ 

In order to calculate -D A1 fcfc(0, 0, 1), we now treat W(x,/i, k) as a small perturbation of W(x, 0,0) for 
|(/it, k)\cx«. "C 1. To this end, we notice that 

8W(n,k)=6W(fi,Q) + k 2 W {fi,k) (18) 

for some matrix Wo(/i, k) which is 0(1) in /i and k 2 . Thus, we must compute 5 W(/x, 0) to first order and 
Wo to leading order in k. The computation of 5W(/j,0) was carried out in detail in [10]. In particular, 
it was proved that if we define the vector solutions corresponding to the functions u x , u a , and ue at 
(//, k) = (0, 0) as Y\, Y 2 , and Y 3 , respectively, then 

/ 

V'(u_)(-T c + u_T a -^T E , 




d_ 

dfi 



n(T, M ,0)l =0 



(9 



and 



^(T,M,0) Uo = 



^3(T, M ,0)U 



du- rT j du— cT , 

— Jo u « da; +wJo 



\ 



V - (l + (c- /'(w-))^-) JoU a dx + (c- /'("-)) %■ uu a dx y 



V - ( 1 + (c - /'0-))^r ) Jo 1 UBda; + ( c ~ W itf uu e<1x I 



where * represents terms which can be explicitly computed but are not necessary at this order in the 
perturbation calculation. Thus, 



d 



SW(/z,0)=<5W(0,0) + Y j (T,fx,0)\^ Q ^e j ] ^ + 0(\fi\ 2 ) 



J =1 



where ej is the j column of the identity matrix on R , and ® denotes the standard tensor product of 
vectors in K 3 . 

Finally, we calculate Wo(0, k) to first order in k 2 . To this end, notice that by the specific forms of the 
first order variations of Y%, I2, and I3 in n, it is enough to compute the k 2 variation in the translation 
direction, i.e. we need only calculate ^p-Yi(T, 0, k) at k = 0. Variation of parameters yields 



^Y 1 (T,0,k)\ k=Q = W(T,0,0) jf WCar.CO)" 1 ^ 





I dx 



du- rT 2 j 
dE JO u x ax 



(c- f'{u-))%- u 2 x dx 
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where the term * is as before. Thus, we see that det (6 W(/x, kj) can be expanded as 

det ^W(0,0)+ ^E^( T '^°)L=o 0e ^ ^+ (^ Y i( T '°' fc )L=o® e i) fc2 j+°(H 4 + fc4 ) ( 19 ) 

in a neighborhood of (/x, k) = (0,0). The matrix involved in the determinant in (TH?)) can be simplified 
via a few elementary row operations. In particular by replacing the third row with the first row times 
the quantity c — f'{uJ) plus the third row, and replacing the first column with the first column minus 
it_/i times the first plus u 2 _fi/2 times the third it follows that the determinant above is equal to that of 
the matrix 

' fc2 W Io U l dx (- TST So Uadx+^ f T UU a dx^j fi (- ^ f Q T U E dx + %jjr UU E dxj fi 

-fN'(uJ)T c V\uJ)T a V'(u_)T E 

\ -fi J Q u a dx -fx J Q u E dx 

A straightforward computation now yields 

i d 3 



det(5W(/i,A)) | (M)=(0;0) - -{T,M} ajE J u 2 x dx 



2 dfidk 2 

as claimed. □ 

Prom Lemma [5] we see that the dominant balance of the equation D(fi, k, 1) — for small \x and k 
is a homogeneous polynomial of degree three in the variables fj, and k. In particular, the three periodic 
eigenvalues at fi — bifurcate analytically from the origin in the transverse wave-number k. If we define 
fi = ak + 0(k 2 ) for |fc| <C 1, it follows that a must be a root of the polynomial 



3 pT 
Q(V) ■= M > P }*,E,c + |{T, M} a ,E J o u\dx. 



Our main theorem on the transverse instability of the periodic traveling wave solutions of the gKdV in 
the ZK equation now follows easily. 

Theorem 1. If {T, M, P} a ,E,c is non-zero at (ao, Eq,cq) € f2, then the spectrally stable periodic traveling 
wave solution u(x;ao, Eq,cq) of ([T|) is spectrally unstable to long wavelength transverse perturbations in 
the gZK equation ^ if {T, M} ajE > at (ao,Eo,co). 

Proof. By solving the equation Q(y) — 0, Lemma [5] implies that there are three periodic eigenvalues in a 
neighborhood of the origin which are given by fiQ = o(k) and 



ilu J ±{T,M} a , E f*u*dx 
^ ±N V {T,M,P} a , E , c + ° (fc) - 

Since u(x;a,E,c) was assumed to be a spectrally stable solution of ((T|), we know from Lemma [T] that 
{T, M, P} a ,E,c > and hence there will be two (non-zero) periodic eigenvalues off the imaginary axis in 
the neighborhood of the origin if {T, M} a ^ E > 0. □ 

Remark 3. Notice that if {T 1 M} aE is negative, Theorem]]] provides no information about the transverse 
instabilities of periodic traveling wave solutions to the gKdV in the gZK equation. This follows from the 
fact that one eigenvalue is just o{k), but we have no information about its confinement to the imaginary 
axis. Moreover, it is clear that Theorem [7] can be extended to the case where {T, M, P} a ,E,c = since 
D(fi,k, 1) is odd in fi and the leading order variation in fi must be negative by Lemma]]] However, the 
vanishing of {T, M, P} a ,E,c will force several branches of spectrum to bifurcate non- analytically from the 
origin. 
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We now point out a few corollaries to Theorem [T] in the case of dnoidal-like solutions for power-law 
nonlinearities f(u) — u p+l , p > 1. By dnoidal-likc solutions we mean solutions which are bounded by 
a homoclinic or heteroclinic orbit in phase space, much like the well known Jacobi dn solutions of the 
focusing mKdV: see [12] for example. To begin, we recall from [10] that the induced scaling in the 
wave-speed yields the following asymptotic relations in a neighborhood of a homoclinic orbitB 

{T, M} a , E ~ -T E M a , {T, M, P} a ^ c ~ -T E M a (— - ±- J P. (20) 



pc 2c 

If we now restrict to dnoidal-like solutions, then, at least for a sufficiently small, it is known by a result 
of Schaaf [22] that T E > (see [10] or [17]). Moreover, it was shown in [10] and [17] that M a < 
for solutions of sufficiently large period, i.e. for orbits sufficiently close to the homoclinic orbit, and \a\ 
sufficiently small. Thus, we immediately have the following corollary of Theorem [T] 

Corollary 1. In the case of a power-law nonlinearity, the dnoidal-like periodic traveling wave solutions 
of (JTJ) of sufficiently large period and \a\ sufficiently small are transversely unstable in the gZK equation 
((2|). Moreover, the periodic traveling wave solutions sufficiently close to a stationary solution in phase 
space are also transversely unstable in the gZK equation. 

Proof. By Lemma [1] and (|20p , such a solution is unstable to long wavelength transverse perturbations in 
the ZK equation if p < 4. Moreover, if p > 4 then {T, M, P} a .E.c < and hence the function D(/i, 0, 1) 
has a non-zero real root /i* > 0. By continuity the function D(fi, k, 1) then, will have a root near /li* 
with positive real part for |fc| <C 1. Thus, although we are not guaranteed a long wavelength transverse 
instability in this case, we still have spectral instability none the less. □ 

In the special case of the KdV equation (/(it) = u 2 ), we are able to use the integrability of the 
traveling wave ODE © to show that {T,M} a . E > for all {a,E,c) <E ft (see [17]). Moreover, up to 
Lie symmetries (Galilean invariance, spatial translation, etc.), all periodic traveling wave solutions of the 
KdV equation can be expressed in terms of Jacobi elliptic functions as 

u(x, t) = 6a 2 7 2 cn 2 (ax + 4a 2 (27 2 - l)t; 7), 

where a > is a scaling parameter and 7 S [0, 1) is the elliptic modulus. In [8], it was shown that all 
such solutions of the KdV are spectrally stable to localized perturbations by using the inverse scattering 
transform. It follows by Lemma[T]that {T, M, P} a ,E,c > in this case. This result can be directly verified 
without the machinery of the inverse scattering transform by using standard elliptic function arguments. 
In particular, noticing that in the case of the KdV the quantities T, M, and P can be represented by 
Abelian integrals of the first, second, and third kinds, respectively, on a Riemann surface it can be shown 
that the Jacobian {T, M, P} a ,E,c can be expressed explicitly as 



rrp M p , _ T 3 (E-V(M/T)) 



where R(u; a, E, c) := E — V(u; a, c) and disc(i?) denotes the discriminant of the polynomial R (see [11] 
for details). Since R(u; a, E, c) clearly has three real roots for each (a, E, c) S ft, and since V(M/T) < E, 
it follows that {T, M, P} a ,E,c > for all (a, E, c) € fi in the case of the KdV equation. Therefore, 
Theorem [l] implies the cnoidal solutions of the KdV are spectrally unstable to long wavelength transverse 
perturbations in the ZK equation. We record this observation in the following corollary. 

Corollary 2. Periodic traveling wave solutions of the KdV with (ao,Eo,co) € SI are spectrally unstable 
to long wavelength transverse perturbations in the ZK equation. 



'Here, we use the notation /(z) ~ g(z) as z — > zq to mean that lim z _> zo f{z)/g(z) = 1 if lim z _> zo g(z) 7^ 0. 
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Finally, we end by analyzing Theorem [T] in the case of the focusing mKdV equation 

u t = u xxx + u 2 u x . (21) 

This equation admits two distinct families of solutions expressible in terms of Jacobi elliptic functions: 
up to Lie symmetries, they are given by 

u{x,t) = V2adn( ) ux + a 2 (2-7 2 );7), (22) 
u(x,t) = V2a~fcn{ux + a 2 (-l + 2j 2 );j). (23) 

In contrast to the case of the KdV, where the elliptic function solutions represented all periodic traveling 
wave solutions, the solutions ([22| and (|23|) only represent the periodic traveling wave solutions of (f2Tj) 
with a = 0. Stability calculations for such solutions has been recently carried out in [11] and [12] in the 
context of nonlinear and spectral stability to various classes of perturbations. Our theory applies not only 
to these elliptic function solutions, but to all periodic traveling wave solutions of the mKdV equation. 
In [11], the signs of the Jacobians {T,M} a ^E and {T, M, P} a ,E,c were computed numerically for all 
(a, E, +1) G 51 using the Abelian integral representation of the necessary quantities. In particular, it was 
shown that for all (a, E, 1) G f2 where {T, M} a ^E is negative, the periodic stability index {T, M, P} a ,E,c 
is also negative and hence one has instability for all (a, E, 1) G fi. We record this observation in the 
following corollary. 

Corollary 3. Periodic traveling waves of the focusing mKdV equation are unstable to transverse pertur- 
bations in the focusing mZK equation. 



4 Transverse Instabilities in the gBBM-ZK equation 

In this section, we consider the spectral instability of a periodic traveling wave solution of the gBBM 
equation to long wavelength transverse perturbations in the gBBM-ZK equation. We begin by briefly 
reviewing the basic properties of the periodic traveling waves of ([3]) (see [18] for more details). For each 
wave speed c > 1, equation (|3|) admits traveling wave solutions whose profiles are solutions of the traveling 
wave ODE 

cu xxt - (c- l)u x + f(u) x = 0, (24) 

i.e. they are stationary solutions of Q in the traveling coordinate frame defined by x — ct. Clearly (|2"4"|) 
defines a Hamiltonian ODE and can be reduced to quadrature: in particular, the traveling wave profiles 
satisfy the relation 

J u 2 + F(u) = au + E 

where a and E are constants of integration and F' = / with F(0) = 0. In order to ensure the existence 
of periodic solutions of (|2"4"1) , we must require that the effective potential 

c — 1 

G(u; a, c) = F(u) — u 2 — au 

have a non-degenerate local minimum. As in the case of the gKdV equation, it follows that the periodic 
solutions of (J2U) comprise a four parameter family of solutions u(x + Xo;a,E,c) while the stationary 
solutions form a codimension two subset. Moreover, we mod out the translation invariance and hence 
only consider the periodic solutions of 1|24|) as consisting of a three parameter family indexed by the 
quantities a, E, and c. 
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The partial differential equation ^ has, in general, three conserved quantities 



T 



M = (u- u xx ) dx 



P = l -j\u 2 + ul)dx (25) 

H = j {^ u2 + F ( u ^j dx 

which correspond to the mass, momentum, and Hamiltonian (energy) of the solution, respectively. These 
three quantities are considered as functions of the traveling wave parameters a, E, and c and their 
gradients with respect to these parameters will play an important role in the foregoing analysis. It is 
important to notice that when restricted to the four-parameter family of periodic traveling wave solutions 
of the mass can be represented as M — u dx. Since all our results concern this four-parameter 
family, we will always work with this simplified expression for the mass functional. 

We now begin our transverse instability analysis in this case. Suppose we have a periodic traveling 
wave solution u(x; a, E, c) of the gBBM equation. We wish to examine the spectral stability of u to long 
wavelength transverse perturbations in the gBBM-ZK equation @. Clearly, u is a y-independent solution 
of the traveling gBBM-ZK equation 

u t - (u xt + u yy ) x - (c - l)u x + f(u) x + cu xxx = 0. (26) 

Linearizing (|26p about u yields the linear partial differential equation 

Vv t = d x (ll[u} + dl) v, 

where T> = 1 — d x and TZ[u] = —cd\ + (c — 1) — f'(u). Seeking separated solutions of the form 

v(x,y,t) = v{x)e^- iky 

where /j, € C and k G R is the transverse wave number leads to the (ODE) spectral problem 

d x (K[u\ + k 2 )v = fWv 

considered on L 2 (R). Notice that the operator T> is positive and invertiblc on L 2 (R). In particular, 
for each fixed k £ R the spectrum of the operator T)~ x d x (lZ[u] + k 2 ) is purely continuous and can be 
described via Floquet theory. In particular, for a given k € R we have that /i £ spec (T>~ 1 d x [TZ[u] + fc 2 )) 
if and only if there exists a A £ S 1 such that 

D{n,k,X) :=det(M(M,fe)- AI) 

vanishes, where M(/i, k) is the monodromy map for the corresponding first order system 

Y x = H(x,fi,k)Y, Y(0, M ,fc)=I, (27) 

where 

( ° 1 

n(x,fi,k) = o o 

V -±(/i + />K) l A c ~ l + k2 -/'(«)) 

We now state the main technical lemma for this section. 

Lemma 3. The equation D(p, k, 1) has the following local normal form in a neighborhood of the origin 
(j*,A) = (0,0): 



.3 



//{T, M, P} a , E ,c + 2fik 2 {T, M} a>E f u\dx + O (|/i| 4 + k q 

Jo 
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Proof. In [18], it was shown that three linearly independent solutions of the equation 7?.[u]w = are given 
by the functions u x , u a , and ue- Letting ~W(x,/j,, k) be a matrix solution of (|27p such that 

cu x cu a CUE 
W(.t,0,0) = CU X x CU ax CUEx 
CU X xx CU aX x CV'Exx 

I T 

and defining S~W(n,k) — W(x, fx, k)\ _., it was shown in [18] that 

D(u,k,l) = --det(5W(fi,k)) . 

c 

and that 

D((m, 0, 1) = —{T, M, P} a , E ,c» 3 + 0(\fi\ 4 ) 

for \/j,\ <C 1. Moreover, the first order variations in fi of the vector solutions Y\, Yz, and Y3 corresponding 
to the functions u x , u a , and ue at /i = were calculated and essentially have the same form as those 
found in the previous section for the gKdV equation. Due to the particular form of i5W(/i,0), it follows 
then that one only needs to compute the first order variation in fc 2 of the vector solution Y\ corresponding 
to the translation direction u x . Using variation of parameters, one can show that 



W(T,0,0)£ W^O.O)" 1 ( 




9u(0) rT 2 j 
C dE Jo u x ax 



l-f'(u(0)))^fj T uldx 



where the term * can be explicitly computed but is not necessary at this order in the perturbation 
argument. A straight forward calculation then gives 

det (fW(ft *)) I (/a)=(0 , ) = - c i T > M }a,E ^ u 2 x dx. 

from which the lemma follows. □ 

Lemma [3] readily yields a necessary condition for the underlying (gBBM) periodic traveling wave 
u(x, a, E, c) to exhibit a modulational transverse instability in the gBBM-ZK equation. 

Theorem 2. If{T, M, P} a ,E,c 7^ 0, then the spectrally stable periodic traveling wave solution of the gBBM 
equation is spectrarlly unstable to long wavelength transverse perturbations in the gBBM-ZK equation if 
{T,M} a , E >0. 

As an application of this theory, one can use the complex analytic methods of Bronski, Johnson, and 
Kapitula [11] in order to prove that in the case of the BBM equation with f(u) = u 2 , one has 

T 3 G' (%) 

{T,M} a , E = — 



12 disc(£- G( • ;a,c)) 



where G is the effective potential energy in this case. In particular, since G' is strictly convex in the case 
of the BBM equation it follows from Jensen's inequality that 



G'(jf\ <±£G'(u(x))dx = 



from which it follows that {T, M} a ^E > for all such solutions. In particular, this proves the following 
corollary of Theorem [2] 

Corollary 4. Letu(x;a,E,c) be a periodic traveling wave solution of the BBM equation. If {T, M, P} a ,E,c 
is non-zero, then u is spectrally unstable to long wavelength transverse perturbations in the BBM-ZK equa- 
tion. 
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Moreover, in [18] the Jacobians {T, M} a ^E and {T, M, P} a ,E,c were analyzed in the case of power-law 
nonlinearities f(u) = u p+l /(p+ 1) in neighborhoods of the homoclinic orbits. In particular, it was shown 
that for all p ^ 4, the quantities {T, M, P} a ,E,c and {T, M} a ^E, c are both positive for all dnoidal type 
solutions of sufficiently long wavelength and \a\ sufficiently small. This prove the following corollary of 
Theorem H 

Corollary 5. In the case of power-law nonlinearities, the dnoidal-like periodic traveling wave solutions 
of ([3]) of sufficiently long wavelength and \a\ sufficiently small are spectrally unstable to long wavelength 
transverse perturbations in the gBBM-ZK equation 

5 Discussion 

In this paper, we extended the methods of [10] to study the transverse instabilities of periodic traveling 
wave solutions of the gKdV equation within a cannonical two-dimensional model: the gZK equation. 
The results are sufficient conditions for the instability to long wavelength transverse perturbations.. In 
particular, it follows from the results in this paper that the cnoidal solutions of KdV are spectrally 
unstable to such perturbations in both the ZK equation, and the dnoidal solutions of the focusing mKdV 
are unstable to such perturbations in the focusing mZK equation. We then extended the results in [18] 
by outlining the corresponding transverse instability theory for the gBBM equation within the gBBM-ZK 
equation yielding analogous results. 

An obvious question arises as to the meaning of the Jacobian {T, M} a ^E in each of these cases. In 
this paper, we have seen that it provides sufficient information for the transverse instability of periodic 
traveling wave solutions of the gKdV and gBBM in several higher dimensional extensions. In [11] and 
[17], this Jacobian was also shown to arise naturally in the orbital stability of such solutions of the 
gKdV to perturbations with the same periodic structure as the underlying periodic wave, and in [18] 
the analogous theory was outlined for the gBBM equation. In particular, in [17] it was found that the 
quantity {T, M} a ^E encodes information about the possible difference in the number of negative periodic 
eigenvalues of the operators C[u] and £[w]|„ , where Hi is the space of mean zero periodic functions: 
notice that by the form of the gKdV all non-trivial evolution of initial data occur within H\. It would 
be interesting to understand the exact physical information the Jacobian {T, M} a ^E encodes about a 
periodic traveling wave solution of gKdV and the gBBM. The key to this may lie in performing a formal 
Whitham theory calculation in the context of slow modulations of a periodic traveling wave solution of 
the gKdV and comparing the formal result to the rigorous Whitham theory type calculatios of this paper 
and [10] and [18]. We have not yet performed this calculation. 

It seems clear that the methods employed in this paper can be applied in a straightforward way to 
studying the transverse instabilities in other nonlinear PDE where the corresponding traveling wave ODE 
is close enough to being completely integrable, in the sense that one can use the integrability to construct 
all elements of the null-space of the linearized operator whether by Noether's theorem or by variation of 
parameters, and where the corresponding linearization yields a spectral problem for a linearized operator. 
It is not clear however if they can be applied to equations such as the gKP equation 

{Ut ~ UXXX ~ f(u) X ) + CTUyy = 

which arises as a weakly two-dimensional variation of the gKdV equation. In particular, the corresponding 
spectral problem arising from linearizing about a periodic traveling wave of the gKdV takes the form 

(dl£[u\ + ak 2 ) v = nd x v 

and hence can not be realized as a spectral problem for any particular linear operator since d x is not 
invertible on Lp er ([0, T]). It would be very interesting to extend the Evans function techniques in this 
paper to consider equations of this form. 
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